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Abstract 

A discussion of the quantum Zeno effect and paradox is given. The quantum Zeno paradox 
claims that a continuously observed system, prepared in a state which is not an eigenstate of 
the Hamiltonian operator, never decays. To recover the classical behavior of unstable systems 
we consider a two-level system interacting with a Bose field, respectively prepared in the excited 
state and in the Poincare invariant vacuum state. Using time-dependent perturbation theory, we 
evaluate for a finite time interval the probability of spontaneous decay of the two-level system. 
Using the standard argument to obtain the quantum Zeno paradox, we consider N measurements 
where iV — > oo and we obtain that the non-decay probability law is a pure exponential, therefore 
recovering the classical behavior. 
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1 Introduction 



In the present paper we are interested to show that a continuously observed quantum system has 
a classical time evolution behavior if it interacts with a vacuum field. Being more specific, we 
study the time evolution of a two-level system, i.e., a qubit, interacting with a Bose field prepared 
in the Poincare invariant vacuum state, showing that the non-decay probability has a exponential 
behavior. 

Unstable systems in quantum mechanics have been the subject of many investigations since 
the origin of this formalism, and have a long story starting with the seminal papers of Gamow [1] 
and Wigner [2]. The temporal evolution of these quantum systems can be roughly divided in three 
different behaviors. A gaussian-like behavior at short times, an exponential decay at intermediate 
times and finally a power-law decay at long times [3] [4] [5]. The gaussian behavior for short times 
is the key point of our discussions bellow, since it leads under general physical conditions to the 
inhibition of the decay of unstable quantum mechanical systems. 

This result was obtained by Misra and Sudarshan [6] [7] using the von Neumann description 
of measurement processes [8]. These authors proved that the realization of many successive mea- 
surements dramatically slow down the evolution of an unstable system. In the limit of continuous 
observations, the temporal evolution of a quantum system can be frozen, defining respectively the 
quantum Zeno effect and paradox, after the paradox given by the Greek philosopher Zeno [9]. 
Actually, many years before, a close related result to the one obtained by Misra and Sudarshan 
was achieved by Khalfin [10], where a proof of the deviation from the exponential decay-law for 
large times was given. 

The first conceptual question in the Misra and Sudarshan construction is the following: we may 
ask whether is possible to realize this limit of continuous observation. Although some authors argue 
that this limit of continuous observation is not physical and must be regarded as a mathematical 
idealization [11] [12] [13], we prefer to discuss this delicate issue later. 

As has been stressed by many authors [14] [15] [16] [17] [18] [19] [20], any quantum mechanical 
system which we are interested to investigate its behavior on time interacts with the surround- 
ings. In other words, in Nature we deal actually with open systems which are influenced by the 
surrounding world through exchange of energy, or in a more abstract way, information. These 
considerations motivate us to investigate unstable systems which are continuously observed in a 
finite time interval, including a vacuum field in the problem. The two-level system, prepared in 
a eigenstate of its free Hamiltonian, can make a transition to the lower eigenstate driven by the 
vacuum field, which acts as a reservoir. This is a development of an old idea. In the derivation 
of Planck's radiation law, Einstein introduced the idea of spontaneous emission, where a system 
makes a transition to a lower eigenstate without external stimulation [21]. 

In quantum mechanics, in treating arbitrary systems the time evolution of closed systems 
is described by one-parameter group of unitary operators, and the equation of motion of such 
systems are symmetric to time reversal. When a small system interacts with an reservoir, which is 
characterized by an infinitely large number of degrees of freedom, the time evolution of the small 
system can not be represented in terms of unitary Hamiltonian dynamics and we call it an open 



1 



quantum system. There are a variety of theoretical models of reservoirs. One situation is when 
the system S is coupled to an infinite number of harmonic oscillators. In this situation there are 
two kinds of reservoir of common interest. The first one is a thermal reservoir, where we assume 
that the harmonic oscillators are in thermal equilibrium at temperature The second one is 
a squeezed reservoir. The specific system-reservoir model which is appropriate for the study of 
several interesting situations is when the harmonic oscillator bath is constituted by a Bose field 
in free space. 

The aim of the paper is to recover the classical time evolution behavior for continuously ob- 
served systems. The situation that we are interested to study is when the two-level system is 
prepared in the excited state and it is interacting with a Bose field prepared in the Poincare in- 
variant vacuum state. Using perturbation theory we compute the probability of decay evaluated 
in a finite time interval. 

These calculations are not new in the literature. Discussing model detectors, Svaiter and 
Svaiter [22] [23] assumed a weak-coupling between a two-level system and a Bose massless field. 
These authors evaluate the transition rates of the two-level system in different kinematic situations 
without use of the rotating- wave approximation [24] [25]. Further, Ford and collaborators [26], 
using the same model, assumed the presence of one or two infinite perfectly reflecting plates 
(mirrors). They show how these mirrors, which change the vacuum fluctuations associated to the 
Bose field, influence radiative processes at zero temperature. They also evaluate the probability 
per unit time of spontaneous emission at finite temperature. Radiative processes of atoms in 
waveguides and cavities also have been investigated by many authors. See for example the Ref. 
[27]. In the Refs. [28] [29] one of the authors continue to investigate radiative processes associated 
to the Unruh-DeWitt detector [30] [31] [32], in interaction with a massless scalar field. Being more 
precise, in Ref. [29], it was calculated the detector's excitation rate when it is uniformly rotating 
around some fixed point, when the scalar field is prepared in the Poicare invariant vacuum state, 
and also when the detector is inertial and the field is prepared in the Trocheries-Takeno vacuum 
state [33] [34]. These two response functions allow to the authors to present questions analogous 
to those discussed by Mach in the Newton's bucket experiment in a quantum mechanical level. 

This paper is organized as follows. In Section II we briefly discuss the theory of the classical 
and quantum mechanical decays and the quantum Zeno paradox. In Section III, assuming that 
the two-level system interact with a Bose field prepared in the Poincare invariant vacuum state, 
the probability of decay of the two-level system is evaluated in a finite time interval. In Section 
IV, using the same arguments as in the quantum Zeno paradox, we obtain that in the case of 
continuous observations the non-decay probability law is exponential for all times. Conclusions 
are given in Section V. In the Appendix we discuss qubit-boson field interaction Hamiltonians. In 
the paper we use ks — c — h — 1 
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2 The classical and quantum mechanical decays and the 



quantum Zeno paradox 



To introduce probability in classical physics we have to make use of a huge number of identical 
prepared systems. The classical theory of decay is quite simple and is based on the assumption 
that unstable systems have a certain probability of decay. The basic features of this simple model 
is that we assume a Markovian approximation. Therefore this probability does not depend on 
the past history of the unstable system. Let us assume N unstable systems, and that the decay 
probability per unit time be a constant that we call T. For simplicity T is characteristic of the 
system and also does not depend on the total number of unstable systems nor on the environment 
surrounding them. Let us define the number of unstable systems at time t by N(t). Therefore the 
number of systems that will decay in the infinitesimal interval of time dt in dN(t). Consequently 
we have 

-dN(t)=NFdt. (1) 

Defining the inverse of T, i.e., the lifetime of the unstable system by te (T = ^), the number of 
unstable systems at a generic time t is 

N(t)=N(0)e X p( r ±), (2) 

where N(0) is the number of unstable systems at the beginning of our observation, i.e., t — 0. 
One define the non-decay classical probability P c iass(t) as 

W*) = g=«p(~> (3) 
For short times (t « t e ) we can write 

Pclass(t) = 1 - — + ... (4) 

Note that we are excluding cooperative effects, therefore Y and also P(t) are environment-independent. 
The solution given by the Eq.(3) has a dissipative behavior and is a fundamental law that gives 
the classical behavior of unstable systems, as in experimental nuclear physics, for instance. 

In quantum mechanics we introduce probability even working with a single system. A quantum 
mechanical treatment of the same problem give us a short and large time behaviors which are in 
disagreement with the exponential law obtained in Eq.(3). Let us first discuss the deviation from 
the exponential decay law for large times. We are following the arguments presented in the Ref. 
[7]. Let us assume a quantum system with a set of observables, i.e., operators which commute 
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with the Hamiltonian of the system. The Hamiltonian has a complete set of eigenstates, which is 
a basis of the Hilbert space, therefore every state vector of the system can be expressed in terms 
of it. For simplicity we assume that there is only one unstable state that we represent by | a), 
which is orthogonal to the bound states of the Hamiltonian operator H . Since the dynamic is 
time-translational invariant, the unitary operator U(t 2 — t x ) propagates the system from t 1 to t 2 . 
Using the self-adjoint operator H of the system, the dynamics is defined by the unitary operator 
U(t) = e~ tHt . Suppose that we are studying the temporal evolution of the system after t = 0. Let 
us define the spectral projection of the Hamiltonian operator by 

H= f dXX\ A) (A|. (5) 

The energy distribution function of the state | a ) or the probability that the energy of the unstable 
state | a ) lies in the interval [E,E + dE] is given by 

r-E+dE 

/ (a\ X) ( A| a) dX. (6) 

JE 

The non-decay probability at the time t is defined by P qua nt(t). Therefore, using the standard 
interpretation of quantum mechanics, we have that 

P q uant(t) = \(a\e- mt I a) | 2 , (7) 

where the decay probability is given [1 — P quan t(t)]- Let us study the non-decay amplitude. It is 
given by 

( a | e~ lHt \a)= f dX e~ iXt ( a\ A) ( A| a) dX. (8) 



If ( A | a) = 0, for A < 0, i.e., the spectrum of the Hamiltonian operator H is bounded from bellow, 
then when t — > oo the quantity P qua ntif) decreases to zero less rapidly than any exponential of the 
form e~ at . Therefore we have a deviation from the exponential decay law at large times. Also, the 
treatment for the same problem for short time give us a short time behavior which is quadratic 
and therefore in disagreement with the exponential law obtained in Eq.(3). 

Let us assume again a quantum system with an Hamiltonian operator H with a complete set 
of eigenstates denoted by ] % ) {% — 1, 2, 3, ...). If we prepare the system in a normalized state | a ) 
which is not an eigenstate of H, it is possible to show that the non-decay probability at short 
times is of the gaussian type. A short time expansion using the Eq.(7) yields 

t 2 

P q uant(t) = 1 - — + (9) 

where the quantity t^ 1 = [( a \ H 2 \a) — ( a | H | a ) 2 ]^ is the inverse of the characteristic time of 
the gaussian evolution. The crucial feature of this approximation is that the non-decay probability 
after a short observation time t is quadratic. The quantity r z is also called the Zeno time. 
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Let us assume that the quantum measurement occurs instantaneously. We also assume that 
is possible to perform infinitely many measurements in a given finite interval. Suppose that we 
perform N measurements at equal time interval which satisfies T = NAt. In each measurement 
we observe that the system stays in the initial state | a ) which was defined before. The probability 
of observing the initial state at the final time T after N measurements reads 



quant 



(T) = 



quant 



N 



Substituting Eq.(9) in Eq.(lO) we have 



quant 



(T) 



1 (T x 2 



For very large N we get 



Pj£L(T) ~ exp 



N 



T 2 



N 



(10) 



(11) 



(12) 



and repeated observations slow-down the evolution of the unstable system and increase the prob- 
ability that the system remains in the initial state at T. If we are able to set N — > oo one 
obtains 

lim P^L{T) » hm exp [-- St 1 = 1. 



(13) 



This is a very simple derivation of the quantum Zeno paradox. The unstable quantum system 
becomes stable if we perform infinitely continuous measurements. 

There are many physical assumptions that we have to make to obtain this effect. Many authors 
claim that the limit of infinite measurements is non-physical, and it is in contradiction with the 
Heiseberg uncertainty principle [35] [36]. We leave open these questions right now and we shall 
come back to this important issue when we discuss the interpretation of time-energy uncertainty 
relations [37] [38] [39] [40] [41]. 

In the next Section we are interested in the quantum measurement of a single object interactng 
with the vacuum modes and how to evaluate the probability of transition in a finite observation 
time. Since we are assuming a weak-coupling between the two- level system and the environment, 
the probability of decay of the two-level system is computed using first-order approximation in 
perturbation theory. 



3 The probability of decay evaluated for a finite time in- 
terval 
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In this section we are interested in computing radiative processes of a quantum two-level 
system, interacting with a vacuum field. For simplicity we will use the following notation. The 
two energy levels of the two- level system, i.e., the ground and excited energy levels, are given by 
uj g and u e (uj e — uo g = uj > 0), with eigenstates of the free two- level system Hamiltonian | g) and 
| e), respectively. We are assuming a non-zero monopole matrix element between these two states 
and we can assume that the diagonal elements of the monopole operator vanish. 

As we discuss in the appendix, the coupling between the massless scalar field and the two-level 
system is given by a monopole interaction Hamiltonian, i.e., 

H I = \m{T)<p{x(r)), (14) 

where m(r) is the monopole operator of the two-level system, (p(x(r)) is the scalar field operator. 
The total hamiltonian of the system is given by Eq.(A.9), the free Hamiltonian of the scalar field 
and the Eq.(A.26), where A is a small coupling constant between the qubit and the quantized Bose 
field. 

We would like to stress that in general, measurement and state preparation are different phe- 
nomena. In quantum mechanics preparing a particular state might involve a special type of 
measurement, but there are an infinite number of prepared states which are not associated with 
measurements. For example, a generic state for the two-level system | can be written as 
| fy) — a | e) + P\ e). The normalization condition gives | a\ 2 + | f3\ 2 = 1. We can also prepare the 
two- level system in another state which is not an eigenstate of the Hamiltonian Hq. Introducing 
the variables 9 and we can write 

| = | 0, 4>) = e _i 2 cos ( ^ J | e) + e*2 sin ( ^ J | g). (15) 
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Clearly \6,4>) is a complete set, since f£ d6sin6 J n ^ | 9, 0) ( 6, <f>\ = 1. As already stated, we 
are interested to study the probability of decay from the excited state, driven by the vacuum 
fluctuations. Therefore we assume that the two-level system is in an eigenstate of the Hamiltonian 
Hq. To evaluate the probability of decay (excitation) of the two-level system interacting with the 
Bose field, we can define the prepared initial state of the system in r = as | r ) = | e ) <g> | <3>j ) 
(I r o ) — I 9 ) ® I ))) where | ) is the initial state of the field. Both situations can be analyzed 
using the same formalism. In the interacting picture, using the first-order approximation the 
probability of transition P(E,t, 0) after the time interval r is given by 

P(E, r, 0) = A 2 | ( e | m(0) | g ) | 2 F(E, r, 0), (16) 

where the response function F(E, r, 0) is given by 

F(E,r,0) = f T dr' f T dT"e- iE ^- T "\^\cp(x(T'))cp(x(T"))\^), (17) 
Jo Jo 
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and in the above equation E — ±cu, where the signs (+) and (— ) represent the excitation and 
decay process, respectively. Note that to obtain Eq.(17) we sum over all possible final states of 
the Bose field and we are using the completeness relation = 1, where | $/) is an 

arbitrary Bose field final state. We are obtaining information about the time evolution of the sub- 
system, i.e., the two-level system. This approach must be equivalent to study the density operator 
of the sub-system, i.e., the two-level system, which describes the dynamic evolution of this sub- 
system interacting with the environment through master equations. The equivalence between the 
density operator description and first-order perturbation theory to evaluate the quantum non- 
decay probability was demonstrated a long time ago by Fonda et al [3] . 

To proceed, let us suppose that it is possible to prepare the scalar field in Poincare invariant 
vacuum state |0,M), or the Minkowski vacuum state. Therefore in the above equation the 
quantity ( <3>i | ip(x(r')) ip(x{r")) \ $j ) becomes the positive Wightman function associated with the 
Bose scalar field evaluated in the world-line of the qubit. There are two points that we would 
like to stress. The first is that in the integrand of Eq.(17), the two-point correlation function 
depends only on the time difference (V — t"). The integration over r' and r" is carried out over 
the square < r' < r, < r' < r [22]. The second one is that we are not using the rotating- 
wave approximation, used by Glauber [42] and others, to define an ideal photo-counter detector. 
Therefore in the response function F(E,r, 0) the vacuum fluctuations contributions associated 
with the Bose field are taken into account and we are studying the radiative processes associated 
to the qubit induced by a vacuum field. 

To calculate the probability of transition evaluated in a finite time interval let us prepare the 
system in the initial instant of time Tj in the state 

\ n ) = \e) ® |0,M), (18) 

and assume that we observe the system in the ground state in the instant of time 77. In this 
situation the response function becomes 

F^(E,T f ,Ti) = r dr 1 P dr" e" 1 ^'-^ (0, M \ <p(x(t')) ip(x(t"))\0, M }, (19) 

J Tj J Ti 

where E = —u. We are using the subscript (1) to call the attention that this quantity is evaluated 
in the first measurement. Defining Ar = 7/ — Tj, and introducing the variables £ = r' — r" and 
i] = t' + t", the response function given by Eq.(19) can be written as: 

f a. (E , AT) = _^£^ (AT _ u|) ^_, (20) 

where the is is introduced to specify correctly the singularities of the Wightman function, to respect 
causality requirements. Let us split the response function F^(E, At) in two contributions: 

F« (E, At) = if } (E, Ar) + F 2 (1) (E, Ar) , (21) 
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where the functions F^\e,At) and F^^E, At) are given respectively by 

^» (B ,Ar) = - i L/^dfAr^J i (22) 

and 

After some calculations [22] [23] we obtain that the functions f[ 1 \e, At) and F^(E, At) can be 
written as: 

We, a,) = £ (-se(-s) + ^ + JfL ( S , £| A, - 0) ,24) 

and 

F 2 (1) (£, Ar) = At (-7 + Ci|£|Ar - In e\E\ - 1) . (25) 
2ir 2 

In the Eq.(24) and Eq.(25), 7 is the Euler constant and the Si(z) and Ci(z) functions are defined 
respectively by [43]: 

Si(z) = [ Z — dt, (26) 
Jo t 

and 

f z 1 

Ci(z) = 7 + lnz+ / -(cosi-l)dt. (27) 
jo t 

The Eq.(25) has two divergences. One given by InAr as Ar — > + and other given by lne. 
In a full perturbative renormalizable quantum field theory, there is a regularization and also a 
renormalization procedure, where the infinities can be eliminated. One way to circumvent this 
problem is to define the rate R^(E,At) = (E 1 , At) [44]. Since we are interested in 

the non-decay probability in a finite time, let us define the renormalized probability of transition 
P£l(E, At) by 

P&(E, Ar) = A 2 |(e | m(r,) | g ) | 2 (f™ (E, At) - In ^) . (28) 

To support this procedure we can use the argument that these divergences are spurious and can not 
appear in the physically measured processes. Using this renormalization procedure, the probability 
of decay can be written as 

P£L(E, Ar) = A 2 | ( e | m(r 4 ) | g ) | 2 F™(E, At), (29) 
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where F r « (E, Ar) is given by 

FW(E, Ar) = JL ^|Ar (| + Si \E\ Ar) + cos £Ar - 1 + j^' ^ 1 (cos£ - 1) d£j . (30) 

For a small time interval Ar the transition probability contains two contributions: the first one 
that increases linearly with the time interval and the second one that increases quadratically with 
the time interval. The probability of non-decay of this two-level system after a finite time interval 
Ar is given by P^ U (E, Ar) = (1 - P«(E, Ar)). 

Instead of using the interaction picture and the perturbation theory in first-order approxima- 
tions, it is possible to use the Heisenberg equations of motion to the Dicke operators and also to 
the annihilation and creation operators associated to the Bose field [45] [46] . Clearly both methods 
of calculations must give identical results. 

In the next Section we use this probability of non-decay to show that the observed probability 
is exponential of the time T. There are some technical problems in the second measurement. 
Although we assume that the interaction between the qubit and the field is weak, we can not 
suppose that the state of the field does not change in time. We conclude that to study the time 
evolution of the system in the second measuremet, we have to assume that the qubit is still in the 
excited state | e), and the Bose field is in an arbitrary state | We will assume for simplicity 
that the initial sate of the field in the second measurement is a many particle state. 



4 The exponential decay after iV successive measurements 

The aim of this section is to show that if we couple the two-level system with the Bose field in 
the vacuum state we recover the exponential non-decay probability if the system is continuously 
observed. 

Suppose that we perform N measurements at equal time interval which satisfies T = iVAr, 
and in each measurement the system stays in the initial excited state | e ) defined before. Again, 
we would like to stress that to obtain the Eq.(17) we assumed that no observation was made 
to discriminate among possible final Bose field states. Therefore to calculate the probability of 
observing the initial state of the qubit | e ) , at the final time T after N measurements we have to 

point out the following fact. In the first measurement, to obtain the probability P^ u (E,At) = 
(l-P£l(E, Ar)), where P£ ] n (E, Ar) is given by Eq.(29) and Eq.(30), we summed over all possible 
final states of the Bose field. This reflect the fact that we are interested in the final state of the 
qubit and not that of the field. For this reason, in the second measurement we can not assume 
that the initial state of the Bose field is identical with the initial state of the field prepared in the 
first time interval Ar. 
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Before the second measurement, the information of the state of the field is retained by the 
system. Since the initial Bose field state in the second measurement is indeterminate for us, we 
choose an arbitrary state | $j ). Therefore to study the time evolution of the system in the second 
measuremet, let us suppose that the qubit is still in the excited state | e), and the Bose field is 
in an arbitrary state | To find the probability of transition from the state | e) ® | $j ) to 
the final state | g ) ® j ) after the second time interval At we have to evaluate the expression 
PW(E, T f + At, n + At) given by 



Note that we are using the same convention used in the section IV. We have again that At = Tf—Ti, 
and £ = r' — r" . Note that to obtain the Eq.(32) we are using again the completeness relation 
over the final states of the field in the second measurement. 

We have a considerable arbitrariness in the choice of the initial state of the field after the first 
measurement. An appropriate starting point is to suppose that the initial state of the field in the 
second measurement is a many particle state with n x quanta with momenta k x and energy ui, n 2 
quanta with momenta k 2 and energy u 2 and so on. Therefore the two-point correlation function 
that appears in the Eq.(32) can be written as 



($, | cp(x(T')) <p(x(t"))\ ) = (n 1 (k 1 )...n i (I 9 )l <P«r')) ^{x{t"))\ n^)..^)). (33) 



Using Eq.(33) it's not difficult to show that we can write the two-point correlation function 
( | (p(x(t')) (p(x(t")) \ $j ) in the following way: 



where G + (x(t'), x{r")) is the positive Wightman function evaluated in the world line of the qubit 
and rii is the number density of quanta in the /c-space. The set (ul_.(x) , u^^x)) is a basis in the 
space of solutions of the Klein-Gordon equation. Without loss of generality we can choose the 
plane-waves for the basis. 

Taking the continuous limit, assuming that the quanta are distributed isotropically, and sub- 
stituting Eq.(33) and Eq.(34) in Eq.(32) we have that the response function F^ 2 \E, At) in the 
second measurement can be written as 




F^(E,At)= [ T d^(AT - \^\)e-^(^\ V (x(T')) V (x(T"))\^}. (32) 



G+(x(t>),x(t")) +J2n i u ki (x(T'))ul(x(T")) +J2n i ul(x(T'))u ki (x(r")), 



(34) 



F i2 \E, At) = Fi 2) {E, At) + F 2 (2) {E, At). 



(35) 
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(2) 

In the Eq.(35), the quantity F{ (E^Ar) is the vacuum contribution given by Eq.(20) and the 
quantity F 2 (E 7 Ar) is the non- vacuum contribution given by 

F?\e,At) = JL j A ^(Ar- ItDe-^g®. (36) 

In the above expression the function g(£) that depends on the the number density of quanta in 
the /c-space in the continuous limit is 

g(£) = J duujn{uj) (e*** + e~ 1 ^) . (37) 

To proceed, we can extend the integration over all frequencies in the Eq.(37) and also replace the 
number density of quanta in the /c-space, n{uj) by a constant value in the interval [0,a]. Using 
that [47] 

L <* ^ ^ (l + £ ^) , 08) 

where / n (x) = x n for < x < a and zero otherwise and n — 1,2, ... we get that Eq.(37) can be 
written as 

g(Z,a) = ~ - 1 (e^ + e-^) - ^ sin£a. (39) 

The first term in the above equation will give a contribution to the response function which is 
proportional to the one of the vacuum field. To proceed, let us substitute Eq.(39) in Eq.(36). 

Therefore we have that F 2 (2) (£, At) can be written as 

F^(E,a,Ar) = jf\E, At) + jf\E, a, At) + fj 2) (E, a, Ar) + ff\E, a, At)), (40) 

where f[ 2 \E, At) = 2F^\e,At), and the two other terms in the above equation are given 
respectively by 

/«<*,«, Ar) = -- 2 (Ar "UD^— (41) 

and fP(E, a, At) = fP(E, -a, At). Again the is is introduced to respect causality requirements. 
Finally, the function /f (E,a, At) is given by 



/i 2, W«.^)^-£ J /_ A >«Ar-UI)e-|^. (42) 
The integral above can be written in the following form: 



ff(E,a,AT) = ( ( AT d A e -^ - f e"^ 



87ri w-Ar £ J-At £ 



a) ) +|^/ o AT ^sin£acos^. (43) 
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To carry out the integrations we use the fact that the two integral in the right side of the Eq.(43) 
can be interpreted as the principal value, i.e. 



e(x) 



-P 



111 



and also the last integral of the Eq.(43) can be carry out immediately and gives 

' (cos(E - a) At - l) — - — - (cos(E + a) At - 1 



I d£sin£acos£P 
Jo 



2(E-a) V 



2(E + a) V 



(44) 



(45) 



We conclude that in this second small time interval At the transition probability contains also two 
contributions: the first one that increases linearly with the time interval and the second one that 
increases quadratically with the time interval. Following this line, the probability of observing the 
initial state at the final time T after N measurements reads 



Using the fact that P^ u (e, ( 



P (1) [ E 

r still 1 ^ 



N 



T. 
N' 



(AT-l) 



(46) 



1-^(^.(5)) 

T 



for i = 1,2, we have 



l-Pi 1 ! E 



N 



l-P£l(E,a, 



T, 
N' 



(JV-l) 



(47) 



For At = j^, expanding for small arguments and keeping terms only through order (J^ the 
quantity F$ n (E, At) is written as 



FW f E, - 

' ' ' ' 1 — _/V 



( \E\T aE 2 T 2 " 
\4ttN + 4nN 2 



(48) 



where a — (| — The quantity F^(E,At) also in the same order of the approximation can 
be written as 



F$ n (E,a, 



T 

N 



p(E,a)?- + q(E,a)?- 



where p(E,a) and q (E, a) are functions of E and a given by 

p(E,a) = ^(3|£| + \E + a\ + \E-a\ +a6(a- E) - aO(-E - a) + aQ[E + a)J 



and 



q(E,a) = ^(5E 2 (7r-3)-2a 2 ). 



(49) 

(50) 
(51) 
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It is easy to see that, if \E\ < \ a\, then 



p(E,a) = ±(3\E\+4a). (52) 

For \E\ > | a\: 

p(E,a) = ^. (53) 

Finally for \E\ = \a\, we obtain 

p(E,a) = ^ (54) 

71 

and 

g(E,a) = (57r-17)^. (55) 

Defining A 2 |( e | m(rj) | g )| 2 = <x, we may write the probability of observation the initial state at a 
finite time T after N measurements P^(E, j^) as 



PsZ\(E,a,N) 



l-^(p(E,a,a)T + q(E,a,a)^r 



(56) 



At this moment we would like to discuss the interpretation of time-uncertainty relations. There 
are a large amount of literature devoted to the interpretation of quantum mechanics. Nevertheless, 
concerning the time-uncertainty relations, there are a few papers discussing the implications of 
such relations. Landau and Peierls [37] and also Landau and Lifshitz [38] claim that the energy of a 
quantum system can be measured exactly at a given time. Nevertheless, we must take into account 
the change caused by the process of measurement. In the relation AE At > 1, the quantity AE 
is the difference between two exactly measured energy values at two different instants of time, 
where At is the time interval between the measurements. If we accept this interpretation, there 
is a finite but very large N constrained by an upper bound given by the Landau, Peierls and 
Lifshitz interpretation of the time-energy uncertainty relation (N < TE), and we obtain that the 
non-decay probability is polinomial but very similar to the exponential behavior. 

On the other hand, this interpretation of the time-energy uncertainty relation is not universally 
accepted. Aharonov and Bohm [40] [41] claim that the time-uncertainty relations are not consistent 
with the general principles of quantum mechanics which require that the uncertainty relations be 
expressible in terms of operators. Therefore, they concluded that the energy of a quantum system 
can be measured in an arbitrary short time. In this framework we are able to take the limit 
N — > oo, and we get 

P^](E,a,T) = exp L p (E,a,a)T-q(E,a,a)^Pj. (57) 
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To conclude, we obtain that repeated observations slow-down the evolution of the unstable system 
and increase the probability that the system remains in the initial state at T. If we are able 
to perform only a finite but very large N number of measurements in a finite time, constrained 
by an upper bound given by the Landau, Peierls and Lifshitz interpretation of the time-energy 
uncertainty relation, we obtain that the non-decay probability is polinomial but very similar to 
the exponential behavior. On the other hand, if we are able to perform continuous observation 
and the limit iV — > oo can be used, the non-decay probability becomes a pure exponential. 

To obtain this classical behavior, we follow different steps. First, we prepare the composed 
system in an initial state | e) <g> | 0, M), i.e., we couple the two-level system with a vacuum field. 
Second, that the system evolves under the influence of the unobserved Bose field. Finally that we 
are able to perform continuous observations in the two-level system. 



In this paper we study the time evolution of unstable systems after repeated but finite observations 
and also in the limit of continuous observation. We show that a continuously observed quantum 
system has a classical time evolution behavior if it interacts with a unobserved vacuum field. 

Using perturbation theory in first-order approximation where a two-level system is interacting 
with the Bose field in the Poincare invariant vacuum state we obtain two distinct types of behaviors. 
First a finite but very large iV constrained by an upper bound given by the Landau, Peierls 
and Lifshitz interpretation of the time-energy uncertainty relation, and second, the case where 
N — > oo, allowed by the Aharanov-Bohm interpretation of the same relation. Studying the non- 
decay probability in both situations, we obtain that the non-decay probability is polynomial and 
very similar to the exponential behavior for the first case. For the second case the non-decay 
probability is given by 



where r c oc Jr. It is important to remind that in the interaction Hamiltonian for our model, we are 
not assuming the rotating-wave-approximation, that excludes terms that represent simultaneous 
qubit and field excitation and de-excitation respectively. Although these terms are not-energy 
conserving, representing virtual processes, we understand that only for very large time interval 
At the contribution coming from these terms can be neglected. In the case that we are interested, 
i.e., a small At, a more carefully procedure is not to omit the energy-non-conserving terms. 
This procedure allow us to obtain the above equation that accounts very well for experimental 
observed facts, as for example the decay of many quantum systems, as unstable atoms or nuclei. 
The result establishes that the quantum theory allow us to recover classical behavior under suitable 
circumstances. 

From the preceding sections it is seen that, although our model of the qubit-Bose field composed 
system is quite satisfactory to obtain the experimentally observed non-decay probability law with 



5 Conclusions 




(58) 
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the exponential behavior, the approach of the paper is intrinsically limited since we would like to 
predict the same behavior in quite general quantum systems. We know that many systems have 
a complete set of discrete eigenstates but also a continuum spectrum. Therefore to construct a 
more realistic model to study radiative processes of unstable systems we have to generalize our 
model to one with two bound states and also a continuum of states. For example in the case 
of the atom, which is a practical photo-detector, there is a continuum of final electron states. 
Assuming the same two-levels and a continuum of states | uj a ), with energy in the range [u c , oo], 
and preparing the small system in the state | e) and again the Bose field in the vacuum state we 
have that the probability, evaluated for a finite time interval, of the system to makes a transition 
to the continuum is given by 



where u ae = u a — u e and p(u a ) is a density of final excited states. Again, the quantity P(ou ae , At) 
is given by 



Note that we have to choose a particular form to the density of final excited states p(u ae ), to 
make sure that the integral given by Eq.(59) converges at infinity. This generalization is under 
investigation by the authors. 

There are also different directions for investigation. To mention a few: first is to assume a 
strong-coupling between the qubit and the Bose field [48] [49]. Second, still in the weak-coupling 
regime, is to assume that the reservoir is in thermal equilibrium or in a squeezed state [50]. It 
can be shown that the behavior of the two-level system in a squeezed bath depends on the way in 
which the squeezed bath is prepared, showing Zeno or anti-Zeno effects. Also it is interesting to 
consider N qubits (N — > oo) interacting with one mode of the Bose field, analyzing the situation 
where the qubit system acts as a reservoir whereas the Bose field is an open system, and study 
the dynamics of the reduced system. 

Finally, this two-level system, referred to as a qubit is the elementary building block of a 
quantum computer [51] [52] [53] [54] [55]. This new area of research has revived the interest 
in open quantum systems. The fundamental technological problem is if it is possible to create 
entanglement properties of states in systems that interact with a reservoir. Several situations 
of entangled systems have been proposed, as for example involving trapping and also cooling a 
small number of atoms. How to isolate atoms from the environment in order to make the effect 
of decoherence negligible is an open problem until now. 

Therefore, another natural extension of this paper is to generalize some results of the paper in 
the case of two- atom systems prepared in an entangled state [56]. Using time-dependent pertur- 
bation theory in a first-order approximation, evaluate the probability per unit-time of decay of the 
symmetric and anti-symmetric states given by Eq.(A.17) and Eq.(A.18) respectively to the ground 
state | gi ) <E> | #2 )■ Note that we have to continue to assume that the qubits also interact with the 
vacuum modes. The possibility to prepare the two-atom system in a entangled decoherence-free 




(59) 




(60) 
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state is a question that has fundamental importance in quantum computing applications. Al- 
though the spontaneous emission from a pair of atoms [57] [58] and the causal aspects of their 
spontaneous decay [59] [60] have been studied by some authors, the extension of this formalism 
for two-atom prepared in an entangled state, evaluating the probability of decay in a finite time 
interval is new in the literature. 
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A The qubit-Bose field interaction Hamiltonians 

In this appendix we consider a very general situation where the system under investigation 
contains a large number of non-identical two-level systems. In order to describe the dynamics 
of the reservoir and the two-level systems we have to introduce the Hamiltonian governing the 
interaction of the quantized Bose field with free qubits. Free means that there is no interaction 
between the qubits. Therefore let us consider a Bose quantum system B, with Hilbert space 
which is coupled with N qubits, with Hilbert space Ti^\ Let us assume that the reservoir is in 
thermal equilibrium at temperature The Bose quantum system is a sub-system of the total 
system living in the tensor product space 

H (B) H (Q) 

Let us denote by H B the Hamiltonian of the quantized Bose field, by Hq the free Hamiltonian 
of the iV-qubits and Hi the Hamiltonian describing the interaction between the quantized Bose 
field and the N qubits. The Hamiltonian for the total system can be written as 

H = H B <g> Iq + Ib <8> Hq + Hj, (A.l) 

where I B and Iq denotes the identities in the Hilbert spaces of the quantized Bose field and the 
N qubits. 

The main purpose of this appendix is to discuss qubit-Bose field interaction Hamiltonians. 
Therefore, let us introduce the Dicke operators to describe each qubit. The free j — th qubit 
Hamiltonian will be denoted by H$, since we are using the Dicke representation. Therefore, we 
have 

where | % }j are orthogonal energy eigenstates accessible to the j — th qubit and uo\^ are the 
respective eigenfrequencies. Using Eq.(A.2) and the orthonormality of the energy eigenstates we 
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can write the j — th qubit Hamiltonian H$ as 



00/ 



(A.3) 



i=l 



Let us define the Dicke operators cr^, cr^ and <r (j) for each qubit by 



and finally 



-0) = § (|2 > < 2j — |1)(1| 
*S) = (12X11),-, 



(A.4) 
(A.5) 

(A.6) 



The Dicke representation is a second quantization of the qubits. Combining Eq.(A.3) and Eq.(A.4), 
the j — th qubit Hamiltonian can be written as 



where the energy gap between the energy eigenstates of the j — th qubit is given by 



, ,(?') , ,(i) 



(A.7) 



(A.8) 



Shifting the zero of energy to \{u)^ + u^) for each qubit, the j — th qubit Hamiltonian given by 
Eq.(A.7) can be rewritten as 

(A.9) 



ffU) — rt(j) n z 



Note that the operators at^, cr,.s and <r^ satisfy the standard angular momentum commutation 
relations corresponding to spin | operators, i.e., 



2a 



_/6 _~r — -T 



and finally 



. <7 (7)' <7 (J) 



-cr 



(A.10) 
(A.ll) 

(A.12) 



A well-known model is a combining system where we have only one mode of the quantized field. 
The Hamiltonian of the j — th qubit Hjp, with the contribution of the one-mode quantized Bose 
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field Hs, and the interaction Hamiltonian can be used to define the Hamiltonian of the total 
system, given by 

I B ® #j? +H B ® I Q + H ( j j) = 

I B <g> a z {j) + uj a^a <g> I Q + g (a + a f ) <g> (cr^ + cr^) , (A.13) 

where the second term in the Eq.(A.13) has the contribution from the quantized Bose field single 
mode Hamiltonian and the last term is the interaction Hamiltonian of the j — th qubit with the 
one-mode quantized field. In the Eq.(A.13) g is a small coupling constant between the qubit and 
the one mode quantized Bose field. The generalization to A" qubits is described by 

N N 

I B ® E H d +h b ® Iq + Y, h \ 3) = 

3=1 3=1 
N N 

3=1 v A j=1 

where the first summation in the right hand side 

E a z {j) = Q w a z {1) ® 1 (8) ... (8) 1 + ... + 1 <g> 1 <8> ... ® 1 (g) fi (JV) trf^, (A.15) 
j=i 

and 1 denotes the identity in the Hilbert space of each qubit. We can also introduce a qubit-qubit 
interaction, which is relevant in the study of entangled states. In an entangled system, the state 
of the composite system can not be factorized in to a product of the states of its sub-systems. For 
example in the case of two-atom systems, takes the form 

2 

In the absence of the " dipole-dipole" interaction the pure Hilbert space of the two qubit system is 
spanned by the states | g± ) <g> | g 2 ), | g\ ) <8> | e 2 ), | e± ) <8> | g 2 ), | e\ ) <S> \ e 2 ), where g and e denotes 
respectively the ground and the excited state. If we include the "dipole-dipole" interaction term 
in the form of Eq.(A.16), the vectors | g\ ) <S> | e 2 ) and | e± ) <S> \g 2 ) are n ot more eigenstates of 
the Hamiltonian of the qubits systems. It can be shown that these two vectors states must be 
substituted by the two entangled states, known in the literature as maximally entangled states 
[56] 

I s ) = -j= (| ei ) ® | g 2 ) + | gi ) ® \ e 2 )) (A.17) 
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and 

\a) = ~j= (|ei) ® \g*) -\gi) ® |e 2 )). (A.18) 

Going back to Eq.(A.14), the interaction Hamiltonian is simplified if we assume the Jaynes- 
Cummings model [61]. Considering the Jaynes-Cummings model for one qubit, we have 

I B ® H ( i ] + H B ® I Q + H\ 3) = 

I B <g> a z (j) + lu a) a ® I Q + g(a® cr^ + a f (g) tr^) . (A.19) 
The generalization to A" qubits is straightforward and is given by 

N N 

3=1 3=1 
N N 

I B ®Y, a (3) + u ° fltfl ® J Q + -7= S ( a ® a J) + flt ® ^O)) • ( A - 20 ) 
j=i viV i=1 

One point which is important to stress is that the terms which we ignore in Eq.(A.19) and Eq.(A.20) 
are the so called counter-rotating terms. This approximation is known as the rotating- wave- 
approximation. In the rotating-wave-approximation we ignore energy non- conserving terms in 
which the emission (absorption) of a quantum of a quantized field is accompanied by the transition 
of one qubit from its lower (upper) to its upper (lower) state. 

So far we have discussed A^ non-identical qubits interacting with one-mode of the quantized 
Bose field. Our aim is now to discuss the interaction of a system of A^ identical qubits with energy 
gap (Q = U2 — UJ1), with an infinite number of harmonic oscillators which defines the reservoir. Let 

b k and bk be the creation and annihilation operators of the k — th harmonic oscillator of frequency 
uok- The total Hamiltonian, i.e., the Hamiltonian of the combined system of the reservoir and the 
N identical qubits interacting with the reservoir reads 

N N 

j=l k vjv j=1 k 

In the Eq.(A.21) the first term is the free Hamiltonian of A^ identical qubits, the second term 
is the free reservoir Hamiltonian and the third term is the interaction Hamiltonian between the 
reservoir and the A^ identical qubits. Notice that we shift the zero of energy for each qubits, as 
we did before, and we are assuming the rotating-wave-approximation, where is the j-th qubit, 
k-th harmonic oscillator coupling constant. 

We can also use a different interaction Hamiltonian, as the one introduced by Di Vicenzo 
[62]. This author proposed a soluble model to study the influence of decoherence in quantum 
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computers, with the following model describing a system of one qubit coupled to a reservoir of 
harmonic oscillators: 

I B <8> H Q + H B <g> I Q + Hj = 

I B <g> Vla z + J2^kb{h ®I Q + gJ2 ( b l + b k) (A.22) 

k k 

where f2 is the usual energy level spacing of the qubit, b\ and b k are respectively the bosonic 
creation and annihilation operators of the harmonic oscillators. Notice the particular coupling 
between the reservoir and the qubit, that allows the loss of quantum coherence induced by the 
reservoir without affecting the qubit. There are two straightforward generalizations for this model. 
The first one is the introduction of a mode- dependent coupling constant [63]. Therefore we have 

I B ® H Q + H B <g> I Q + Hi = 

I B ® ft o z + b{b k (g) I Q + ( X k b\ + K h) ® ° z ■ (A.23) 

k k 

Other straightforward generalization is to introduce identical qubits and the Hamiltonian of 
the composed system reads 

I B <g) Hq + H B Iq + Hi — 




(A.24) 



Another generalization of the Hamiltonian given by Eq.(A.21) is not assume the rotating- wave- 
approximation in the interaction Hamiltonian. Without the rotating-wave-approximation, the 
interaction Hamiltonian between the qubits and the reservoir of harmonic oscillators reads 

N 

Hi = -j= E E + b l) ® (-5) + %)) ■ {A.25) 

To conclude this appendix we should point out that through the paper we used the simple model 
where the interaction Hamiltonian between the two-level system and the scalar field is linear in 
both field and qubit and is given by 

Hi = X (m 2 i a + + m 12 a" + a z (m 2 2 - m u )) <g> ip(x), (A.26) 

where — ( i \ m(0) \j), and A is a small coupling constant. The Bose field (p(x) can be expanded 

as 

<p(?) = E ( a *Mt, x) + 4<(*. x )) • ( A - 27 ) 

k 
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The modes Wk(£,x) form a basis in the space of solutions of the Klein-Gordon equation. It's 
convenient to restrict the u^(t, x) to the interior of a three dimensional torus of side L (i.e., choose 
periodic boundary conditions). Then 

u k (t,x) = (2L 3 cu)- 1/2 e ik - x -^ (A.28) 

where 

ki = 2nji/L ji = 0, ±1, ±2... i = 1,2, 3. (A.29) 

It is possible to show that the interaction Hamiltonian defined by Eq.(A.26) is equivalent to the 
interaction hamiltonian given by Hj = m(r)</?(x(r)). This model is known as the Unruh-Dewitt 
detector. The detector is an idealized point-like object with internal degrees of freedom defining 
two energy levels. Different coupling between the field and the two-level system was analyzed by 
Hinton [64]. 
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